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Abstract 

We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic simulations of the 
rotational collapse of a supernova core. In order to calculate the waveforms of the gravitational 
wave, we derive the quadrupole formula including the contributions from the electromagnetic fields. 
Recent stellar evolution calculations imply that the magnetic fields of the toroidal components are 
much stronger than those of the poloidal ones at the presupernova stage. Thus, we systematically 
investigate the effects of the toroidal magnetic fields on the amplitudes and waveforms of the 
gravitational wave. Furthermore, we employ the two kinds of the realistic equation of states, which 
are often used in the supernova simulations. Then, we investigate the effects of the equation of 
states on the gravitational wave signals. As for the microphysics, we took into account electron 
capture and neutrino transport by the so-called leakage scheme. With these computations, we 
find that the peak amplitudes of the gravitational wave are lowered by an order of 10% for the 
models with the strongest toroidal magnetic fields. However, the peak amplitudes are mostly 
within the sensitivity range of laser interferometers such as TAMA and the first LIGO for a source 
at a distance of 10 kpc. Furthermore, we point out that the amplitudes of second peaks are still 
within the detection limit of the first LIGO for the source, although the characteristics of second 
peaks are reduced by the magnetic fields. We stress the importance of the detection, since it will 
give us information about the angular momentum distribution of massive evolved stars. When we 
compare the gravitational waves from the two realistic equation of states, significant differences 
are not found, except that the typical frequencies of the gravitational wave become slightly higher 
for the softer equation of state. 

PACS numbers: 04.30.Db 
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I. INTRODUCTION 



Rotation has been supposed to play an important role in the gravitational radiation from 
core collapse supernovae. The large-scale asphericities at core bounce induced by rotation 
can convert the part of the gravitational energy into the form of the gravitational waves. 
The expected amplitude of the gravitational waves from a supernova in the Milky Way 
is considered to be within the detection limits of the long-base line laser interferometers 
such as [GEO600, LIGO, TAMA, VIRGO] The detection of the gravitational wave 
is important not only for itself but also for the understanding of core collapse supernovae 
themselves, because the gravitational wave is the only window that enables us to see directly 
the innermost part of an evolved massive stars, in which the angular momentum distribution 
and the equation of state are uncertain. 

So far there has been extensive work devoted to studying gravitational radiation in ro- 
ta«o„a, CO. coUap. S fl B B B B 3 (.ee Q fo. a ™, Mo. .cenU. O. 
al. jl^ performed a large number of purely-hydrodynamic calculations employing a realistic 
equation of state (EOS) while neutrino transfer and microphysics are not treated. They 
investigated the effects of initial rotation rates and degree of differential rotation, on the 
g™v,tat,o.al wavefo™. (see also flB). 0„ the othe. hand. MUlle. et aL Q perfo^ed a 
small set of the rotational core collapse simulations while employing the elaborate neutrino 



transport with the detailed microphysics. Here it is noted that these recent [11|, studies 
took into account the initial models based on the recent stellar evolution calculations jl^j ]. 
which predict rather slow rotation rates at a presupernova stage than those assumed in the 
previous studies. If so, the peak amplitudes may not be obtained at core bounce but at the 
ater phases when the neutrino-driven convections become active behind the accreting shock 



121. 



In addition to rotation, we study the effects of magnetic fields on the gravitational signals 
in this paper. In order to estimate the gravitational waveforms in the magnetohydrodynamic 
computations, we modify the conventional quadrupole formula including contributions from 
the electromagnetic fields. To our knowledge, the effects have not been investigated so 
far. To be rigorous, it is true that the realistic magnetorotational core collapse simulations 
should require the implementation of a realistic EOS, an adequate treatment of microphysics 
(weak interactions with neutrino transfer), and relativistic treatment of gravitation in three 
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dimensions. However it is far beyond our scope to treat them all at once. Thus, we choose to 
employ a realistic EOS and treat the microphysics in a simplified manner in the Newtonian 
gravity assuming axial symmetry. Recent stellar evolution studies show that the toroidal 
ma gnet ic field components may be stronger than the poloidal ones prior to the collapse 
[islll^. We systematically change the strength of rotation and toroidal magnetic fields in 
a parametric way. By so doing, we hope to understand the effect of toroidal magnetic fields 
on the waveforms and the amplitudes of the gravitational wave both in the weak magnetic 
field (~ 10^^ G) regime and in the strong magnetic field (~ 10^^ G) regime at the formation 
of protoneutron star. The latter case may be associated with the formation of the so-called 
such as anomalous X-ray pulsars and soft gamma-ray repeaters [16|, 



magnet ars 



Although the number of the magnetars is much smaller in comparison with the canonical 
pulsars, the gravitational wave from such objects should be investigated. 

Furthermore, we investigate the effects of realistic EOS's on the gravitational signals. 
Needless to say, EOS is an important microphysical ingredient for determining the dynamics 
of core collapse, eventually, the gravitational wave amplitude. As a realistic EOS, Lattimer- 
Swesty EOS [\M has been used in recent simulations discussing gravitational radiation from 

nn 

rotational core collapse I12|. It has been difficult to investigate the effects of the EOS's 
on the gravitational si gna ls because available EOS's based on different nuclear models are 



limited (see, however, 
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ulations has become available 
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). Recently, a new complete EOS for supernova sim- 



The EOS is based on the relativistic mean field 
(RMF) theory with a Thomas- Fermi approach. By implementing these realistic EOS's into 
the magnetohydrodynamic simulations, we first investigate the effects of the realistic EOS's 
on the gravitational wave signals. 

We describe the numerical models in the next section. In the third section, we show the 
numerical results. A conclusion is given in the last section. 

II. MODELS AND NUMERICAL METHODS 
A. Initial Models 

Recently, Heger et al. performed the stellar evolution calculations , in which rotation 
and magnetic fields are taken into account. Remembering caveats that their calculations are 
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based on the one- dimensional models with uncertainties and not the final answer probably, 
they pointed out that the toroidal components of the magnetic fields are stronger than the 
poloidal ones at the presupernova stage. It is because the differential rotation amplifies the 
toroidal components by a convective stellar dynamo during the quasistatic stellar evolution 



|26[. To our knowledge, most of the past MHD simulations for investigating the 
dynamics of core collapse supernovae chosen the poloidal magnetic fields as the initial con- 
ditions. This situation motivates us to take the initial conditions, in which the toroidal 
c„,npone.ts a.e do,„,.a„t ove. tKe po.o.dal o.es. S„.e Hege. . al Q calculated o,U. a 
small set of models so far, we prefer a parametric approach to construct the initial conditions 
in this paper. 

We assume the following two rotation laws. In addition, we prepare the toroidal magnetic 
fields to yield the same profile as the rotation. 
1. Shell-type rotation: 
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m = fio X (1) 

Mr) = Bo X (2) 

where Q{r) is angular velocity, B^{r) is toroidal component of the magnetic fields, r is radius, 
and Qq, Rq are model constants, 
2. Cylindrical rotation: 

n(X.Z)^noX3^.^. (3) 

i5*(A-.Z)^B„x^^.^. (4) 

where X and Z denote distances from the rotational axis and the equatorial plane, and 
Xq, Zo are model constants. The other parameters have the same meanings as above. 

We have computed 14 models changing the combination of the total angular momentum, 
the rotation law, the degree of differential rotation, the total magnetic energy, and the 
equation of state. The model parameters are presented in Table H] The models are named 
after this combination, with the first letter, S (slow), M (moderate), R (rapid) representing 
the initial T/|iy|, the second letter S (shell-type), C (cylindrical), denoting the rotation 
laws, the third letter, L (Long), S (Short), indicating the values of Rq, Zq, which represent 
the degree of differential rotation, and the fourth letter, 7 ~ 0.3, indicating the value of 
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-E-m/lW^I- It is noted that the ratio of magnetic and rotational energies to gravitational 
energy are designated as i?m/|W^| and T/|iy|, respectively. We have chosen ~ 10"'^, 0.5, 1.5% 
for the initial Tl\W\ and ~ 10"^ 10-^ IQ-^ lO'S 10-°-3(~ 0.5) % for the initial E^/\W\. 
The initial poloidal components of magnetic fields are assumed to be uniform and parallel 
to the rotation axis, whose value is taken to be about four orders of magnitude lower than 
the toroidal components in accordance with the results by Heger et al. 'l3|. We employ the 
Lattimer-Swesty EOS in Model MSL4-LS, on the other hand, the relativistic EOS in the 
rest of the models. We have made precollapse models by taking a density, internal energy, 
electron fraction distributions from the spherically symmetric 15Mq model by Woosley and 
Weaver 2^ and adding the angular momentum and the magnetic field according to the 
prescription stated above. 

Heger et al. pointed out that the iron core may rotate more slowly with the magnetic 
fields than without. It is because the magnetic braking reduces the angular momentum of 
the core during the quasistatic stellar evolution. Model SSL7 corresponds to the magnetoro- 
tational progenitor by them. We note that the models with the strongest magnetic fields are 
probably unrealistic as suggested from the results by Heger et al. [13]. However we prepared 
these models in order to cover the wide range of the field strength and to see the effects of 
the magnetic fields on the gravitational wave signals clearly. 



B. Magnetohydrodynamics 

The numerical method for MHD computations employed in this paper is based on the 
ZEUS-2D code j2^. The basic evolution equations are written as follows, 

| + PV..^0. (5) 

= -VP-pV$ + ^(V X B) X B, (6) 



P—JT- = -^V ■ V, (7) 
V x{vxB) (8) 



dt 
OB 



dt 

A$ = AttGp, (9) 
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TABLE I: The model parameters. 



Model 


T/\W\{%) E,^/\W\{%) Rotation Law Rq, Xq, Zq xlO^ 


(cm 


-) ^^0 (s-^) 


Bo (G) 


SSL7 


3.2 X 10-3 


9.7 X 10-^ 


Shell-type 




Rq = 1 




0.1 


5.0 X 10^ 


MSL4 


5.0 X 10-1 


10-4 


Shell-type 




Ro = l 




4.0 


1.6 X lOii 


MSL2 


5.0 X 10-1 


10-2 


Shell-type 




Ro = 1 




4.0 


1.6 X 10i2 


MSLl 


5.0 X 10-1 


10-1 


Shell-type 




Ro = 1 




4.0 


5.0 X 1012 


MSS4 


5.0 X 10-1 


10-4 


Shell-type 




Ro = 0.1 




63.4 


3.5 X 1012 


MSS2 


5.0 X 10-1 


10-2 


Shell-type 




Ro = 0.1 




63.4 


3.5 X 10i3 


MSSl 


5.0 X 10-1 


10-1 


Shell-type 




Ro = 0.1 




63.4 


1.1 X 10i4 


MCS4 


5.0 X 10-1 


10-4 


Cylindrical 


^0 


= 0.1,Zo = 


1 


44.4 


1.5 X 1012 


MCS2 


5.0 X 10-1 


10-2 


Cylindrical 


^0 


= 0.1,Zo = 


1 


44.4 


1.5 X 1013 


MCSl 


5.0 X 10-1 


10-1 


Cylindrical 


^0 


= 0.1,Zo = 


1 


44.4 


4.8 X 10i3 


MCS0.3 


5.0 X 10-1 


10-0-3 ~ 0.5 


Cylindrical 


^0 


= 0.1,Zo = 


1 


44.4 


1.0 X 1014 


RCSl 


1.5 


10-1 


Cylindrical 


^0 


= 0.1,Zo = 


1 


76.8 


4.8 X 1013 


RCS0.3 


1.5 


10-0-3 ~ 0.5 


Cylindrical 


^0 


= 0.1,Zo = 


1 


76.8 


1.0 X 10i4 


MSL4-LS 5.0 X lO'^ 


10-4 


Shell-type 




Ro = l 




4.0 


1.6 X lOii 



where p, v, e, P, B, $ are density, velocity, internal energy, pressure, magnetic field, gravita- 
tional potential, respectively. We denote the Lagrangian derivative as d/dt . The ZEUS-2D 
is an Eulerian code based on the finite-difference method and employs an artificial viscosity 
of von Neumann and Richtmyer to capture shocks. The time evolution of magnetic field 
is solved by the induction equation, Eq. (jHI). In so doing, the code utilizes the so-called 
constrained transport (CT) method, which ensures the divergence free (V • -B = 0) of the 
numerically evolved magnetic fields at all times. Furthermore, the method of characteristics 
(MOC) is implemented to propagate accurately all modes of MHD waves. The self-gravity 
is managed by solving the Poisson equation with the incomplete Cholesky decomposition 
conjugate gradient method. In all the computations, spherical coordinates are used and 
one quadrant of the meridian section is covered with 300 (r) x 50 {$) mesh points. We 
made several major changes to the base code to include the microphysics. First, we added 
an equation for the electron fraction to treat electron captures and neutrino transport by 
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the so-called leakage scheme 129 
separately from the main hydroc 
implemented a relativistic EOS 



30 



31 



22] • The calculation of electron fraction is done 



mamic step in an operator-splitting manner. Second, we 



or Lattimer-Swesty EOS [3] instead of the ideal gas 
EOS assumed in the original code. For a more detailed description of the methods, see 
Kotake et al. js^. 



C. Gravitational wave signal from the magnetized stellar cores 

Applying the methods 3,0,131, we derive the quadrupole formula in order to compute the 
gravitational waveforms from the magnetized stellar cores. The dimensionless gravitational 
wave amplitude h^^ = g^^ — rj^y can be calculated by the quadrupole formula as follows: 

(10) 

where i,j run from 1 to 3, t is the time, R is the distance from the source to the observer, 
the superscript "TT" means to take the transverse traceless part, and Uj is the reduced 
quadrupole defined as 

hj = J p*{x,t)^XiXj - ^x^Sij^ d^x, (11) 

where p* represents the total energy density including the contribution from the magnetic 
field, 

P. = P+gf^. (12) 

Expanding the right hand side of Eq. (|lUj) in terms of the pure-spin tensor harmonics 
assuming axial symmetry shows that there is one nonvanishing quadrupole term, namely 
A2Q. Then one derives for the component of h'^^ the following formula, 

where a is the angle between the symmetry axis and the line of the sight of the observer. 
Here Afg is defined by the second time derivative of the mass quadrupole formula: 

= ^M,-o^ (14) 

where the mass quadrupole formula is given as 

G'327r3/2 1 



E2 

20 - ^4 



f°° /3 1\ 



where n = cos 6. By a straightforward, however tedious, calculation to replace the time 
derivatives by the spatial derivatives applying the continuity equation, Eq. (0), the equation 
of motion, Eq. ©, and the induction equation, Eq. (jSI) noting the divergence- free constraint 
(V ■ -B = 0), can be transformed into the following form. 



aE2 ^ aE2 . aE2 

^20 — ^20 ,quad ' ^20 ,Mag' 



where Afg ^^^^ is the contribution from the matter: 



(16) 



aE2 

^20 ,quad 



G 327r3/2 / 1 



dfl J r"^ dr p[Vr'^{3f/ — 1) + Vg'^{2 — 3fl^) — — 6VrVg fl\jl — IJ? 
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-ra^$(3/i2 - 1) + 390$ /iJl - ;u2] + 



djj, 1 r dr\qr{3j^ — 1) — 3 qe fiJl — jJ? 
Jo ^ 



(17) 



AE2 

^20 Mag 



= A 



E2 

20 jxB 



+ p is the contribution from the magnetic field: 



E2 _ G'327rV2 1 



r (ir 



(3/1^ - 1) - (j X S), - 

c 



3^^1 - /i2 - (j X S) 



(18) 



^E2 
^20 p„ 



(iyU / dr 

15 "^0 JO Sucdt 



^[By{3 fi^ - 1)]E^ - —[Ber\3fi^ - l)]rE^ + 



+ ^[B^r'{3f,'-l)]rEg-^^ 
or sm 9 o9 



[B^sm9r%3fi'^ - l)]Er 



(19) 



^20 jxB' ^20p represent the contribution from j x B part and from the time derivatives of 
the energy density of electro-magnetic fields, respectively. We take the first time derivative 
of the magnetic fields, because this is the leading order and the numerical treatments of 
the second time derivatives are formidable. Here, we write the gravitational amplitude as 
follows for later convenience, 



h 



TT 



"-quad ^ "'jxB ^ 



Pm ' 



(20) 



where the quantities of the right hand of the equation are defined by Eqs. (IT!?|) . (fT7j) . (|TH|l . 
and (fT^ . It is noted that we take into account the terms related to the artificial viscosity 
(see Eq. (fTTj) . e.g., 2\). As mentioned, we employ an artificial viscosity of von Neumann 
and Richtmyer. The concrete form of is. 



(21) 



where i = r,6 with / defining the dissipation length. In the following computations, we 
assume that observer is located in the equatorial plane since the gravitational wave radiates 
most in the plane (a = n/2 in Eq. (fT^ ). and that the source is assumed to be located at 
our galactic center {R = 10 kpc). 



III. NUMERICAL RESULTS 

A. Effects of the magnetic fields 

We will first show the effect of the magnetic fields on the amplitude of the gravitational 
wave. For later convenience, the values of several important quantities are summarized in 
Tables ITTl and UTTl We find that the amplitude is affected in the strongly magnetized models 
whose initial i?m/|iy| is greater than 0.1 %. It is natural because the amplitude contributed 
from the magnetic fields should be an order of -Rmag in comparison with the mass quadrupole 
moment component (see Eq. ()12|l ). where 

Bl B 2 

i?mag = ^ ~ 10 % (several x lO^^ q) (lO^^ g'cm^s) ' (22) 

with Be, pc being the central magnetic field and the central density. Thus, strongly mag- 
netized models, whose central magnetic fields at core bounce become as high as ~ 10^'' G, 
can affect the amplitude. As for the waveforms, we find that the contribution from j x B 
part dominates over one from the time derivatives of the energy density of electro-magnetic 
fields (see the left panel of Figure ^). Furthermore, it is found that the j x B part changes 
at the opposite phase of the matter quadrupole components (see the right panel of Figure 
HJ. As a result, the negative parts of the amplitudes become less negative, while the positive 
parts become more positive. This lowers the peak amplitude at core bounce by an order of 
10 % (see Table HTH) . If the initial strength of the magnetic field is the same, the effect of 
the magnetic fields on the amplitude becomes more prominent for the fast rotation models 
(compare the values of ratio of MCSl with RCSl, and MCS0.3 with RCS0.3 in Table ITTH) . 
This is because the central density at core bounce is lowered by the faster rotation (see Table 
HHandEq. 

For the models, whose initial E'm/lVTl's are less than 0.1%, the maximum amplitude is 
found to be determined by rotation as seen from Table |n] We compare the waveforms for the 
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models MSL4 (left panel), MCS4 (right panel) in Figure |21 If we follow the commonly used 
categories of the waveforms 0, |^ , the left panel shows type I behavior (relatively shorter 
intervals of the spikes) and the right panel type II behavior (longer intervals of the spikes). 
What determines the difference is the initial rotation rate and degree of differential rotation. 
It is found that type II occurs for the models with the strong differential rotation with a 

n 

cylindrical rotation law as in the pure rotation case by Kotake et al. |35:] and Ott et al. 



Furthermore, the sign of the values of the second peaks is found to be negative for the 



strongly differentially rotating models with a cylindrical rotation law, and positive for the 
other models (see Table HVjl . Note that by the "second peak" we mean where the absolute 
amplitude is second largest. Due to the j x i? part, the second peaks with the positive values 
become more positive, while those with the negative values become generally less negative 
(see Table CVl). 



B. Effects of the equation of states 

We compare the models MSL4 with MSL4-LS in this section. We repeat that these 
models differ from its employed EOS's. Other than this, any conditions, such as rotation, 
magnetic field, and the employed microphysics, are the same. Therefore the differences 
in the following purely reflect the influences of the employed EOS's. The most important 
difference of the two EOS's is the stiffness of the EOS's at nuclear matter density. In fact, 
the value of the incompressibility of the relativistic EOS (K = 281 MeV) is larger than 
that of Lattimer-Swesty EOS (K = 180 MeV). This means that the Lattimer-Swesty EOS is 
softer than the relativistic EOS. We note that there are three choices with different values of 
(K = 180, 220, 375 MeV) in the Lattimer-Swesty EOS. We take the most softest one in this 
calculation. This is because the value has been often employed for core collapse supernova 
studies. 

A core bounce occurs when the central density reaches its peak of 3.2 x lO^'* g cm~^, 
2.6 X 10^^ g cm-3 at = 215 ms, 243 ms for the models MSL4-LS, MSL4, respectively. 
The earlier core bounce with the higher central density for the model MSL4-LS is due to 
the softness of the EOS. The softness of the Lattimer-Swesty EOS can be seen from Figure 
ini The soft EOS results in the more smaller inner core of Mic = 0.69Mq than that of 
Mic = O.83M0 by the relativistic EOS (see Table ITT|) . This result can be understood as 



11 



follows. The mass of inner core at core bounce is proportional to the square of the lepton 
fraction, YJ. The electron capture with the neutrino emission proceeds further for the soft 
EOS because the core can contract more compact. As a result, the lepton fraction at core 
bounce becomes smaller, which results in the smaller mass of the inner core. 

In the left panel of Figure ^ the gravitational waveforms for the models are presented. 
Shorter burst intervals for the model MSL4-LS are also due to the softness of the EOS. 
The duration of the burst are related to the typical dynamical timescale Tdyn ~ 1/ \fGp^i 
where pc is the central density. Since the values of pc at the core bounce and the subsequent 
oscillations are larger for the soft EOS (see the right panel of Figure |^, the intervals of the 
burst become short. Although the collapse dynamics is changed by the differences of the 
EOS's, it is found that the values of the maximum amplitudes for the two models do not 
differ significantly (see Table ITT|l . This can be understood as follows. The amplitude of a 
gravitational wave is roughly proportional to the inverse square of the typical dynamical time 
scale, (see Eq. (|TO|l ). Since tdyn is proportional to the inverse root of the central density, the 
amplitude becomes larger for the soft EOS by this factor. On the other hand, the amplitude 
is proportional to the value of the quadrupole moment, which becomes in turn small for 
the soft EOS due to the smaller inner core. The maximum amplitude is determined by the 
competition between these factors. As a result, the amplitude becomes almost the same for 
the Lattimer-Swesty EOS and the relativistic EOS. 

In our calculations, we note that there are no models that correspond to the so-called 
type III waveform, which occurs only when a EOS is very soft before core bounce jj]. This 
is because the relativistic EOS and Lattimer-Swesty EOS is not so soft in the corresponding 
regime. 



C. The gravitational wave properties in the weakly magnetized and slowly rotating 
model 

We show the properties of the waveform with collapse dynamics in the model SSL7. It 
is noted that the initial condition for the model is based on the recent stellar evolution 
calculation [3]. We repeat that the initial and -Eni/|W^| are much smaller than the 

other models. The model bounces at t = 215 msec at a central density of 3.0 x lO^^g cm~'^. It 
occurs not by rotation or magnetic fields but by the stiffening of the equation of state. After 
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the core bounce, very weak signals are lasting (see Fig. EI). The gravitational amplitudes 
from the bounce and later reexpansion phases are three orders of magnitude smaller than 
the other models. The timescale of the spikes (~ msec) is determined by the small scale 
non-radial motions behind the shock wave. 

D. Maximum amplitude and second peak 

In Figure IHl the absolute values of the maximum amplitudes and second peaks are pre- 
sented. Note that the data points are confined to small regions although we explore the 
wide range of the initial magnetic field strength. This is because the magnetic fields lower 
the peak ampitudes by ~ 10% but do not change the typical frequencies of the gravitational 
waves significantly. Thus the maximum amplitudes and the typical frequencies are mainly 
determined by rotation. The maximum amplitudes are mainly clustered in the two regions 
(see the open squares in Figure E]). One is around ~ 100 Hz for the models whose names 
begin with "R" (rapidly rotating models) and another is around from ~ 300 Hz to ~ 500 Hz 
for the models whose names begin with "M" (moderately rotating models). The values of 
the maximum amplitudes are in the range of 3.5 x 10^^^ < h'^^ < 2.6 x 10"^'^. The smallest 
value is from model SSL7, whose initial model is taken from the recent stellar evolution cal- 
culation. If the model is correct, it should be difficult to detect the gravitational wave even 
for the next generation detectors such as advanced LIGO and LCGT unless a supernova 
occurs very close to us (see Figure IHl)- When we compare the maximum amplitudes from 
model MSL4 (relativistic EOS) with that from model MSL4-LS (Lattimer-Swesty EOS), 
no significant differences are found, except that the typical frequency for LS EOS becomes 
slightly higher due to the softness of the EOS (see Figure Ej). Therefore the features of the 
gravitational wave are almost independent of the realistic EOS's. 

As pointed out earlier, we found that the signs of the values of the second peaks are 
negative for models with strongly differential rotation with a cylindrical rotation law and 
positive for the others (see Table IIII|) . The absolute amplitudes of the second peak are 
also presented in Figure El As shown, they are within the detection limit of first LIGO for 
a source at a distance of 10 kpc, although the absolute values of the negative values are 
reduced by the magnetic fields. In addition, it is quite likely that the detectors in the next 
generation, such as advanced LIGO and LCGT, to detect the difference. We will obtain 
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information about the angular momentum distribution of evolved massive stars, if we can 
detect the differences by observations of the gravitational wave. However, it seems difficult 
to get information about the distributions of the magnetic fields only by the gravitational 
wave from the magnetorotational core collapse. We note that the effects of the magnetic 
fields on the signals rnay become important after the formation of the strongly magnetized 



star. As pointed 



40 



41[, the gravitational waves, which are detectable for interferometers 



such as LIGO, may be emitted by the global rearrangement of the strongly magnetic fields 
{B > lO^^G) in the core. 



IV. CONCLUSION 



We have done a series of two-dimensional magnetohydrodynamic simulations of the 
rotational collapse of a supernova core and calculated gravitational waveforms. We have 
modified the conventional quadrupole formula in order to include the contributions from 
the electromagnetic fields. Recent stellar evolution calculations imply that the magnetic 
field of the toroidal component is much stronger than that of the poloidal ones at the 
presupernova stage. In this study, we systematically investigate the effects of the toroidal 
magnetic fields on the gravitational waveforms. We have employed the two kinds of the 
realistic equation of states, which are often used in the supernova simulations. By so doing, 
we have investigated the effects of the equation of sates on the gravitational signals. As 
for the microphysics, we took into account electron capture and neutrino transport in an 
approximate method. We found the following. 

(1) . Effects of the magnetic field on the amplitude and the waveform of the gravitational 
waves appear for the strongly magnetized models, whose initial i?m/|iy|'s are greater than 
0.1%. As for the contributions from the magnetic field, the j x B part is found to dominate 
over the time derivatives of the magnetic energy. Since j x B part changes with the op- 
posite phase in comparison with the mass quadrupole moment components, the maximum 
amplitudes are lowered by an order of 10% for the strongly magnetized models. However, 
the maximum amplitudes are mostly within the detection limits of the detectors of TAMA 
and first LIGO if a source is located at a distance of 10 kpc. 

(2) The maximum amplitudes do not change significantly if we employ the two realistic 
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TABLE II: Summary of important quantities for all models, tb is the time of bounce, p^axh is the 
maximum density at bounce, Mj.c b is the mass of inner core at bounce, r/|VF|finab -f'm/lW^lfinai is 
the final ratio of rotational and magnetic energies to gravitational energy of the core, respectively. 
At is the duration time full width at half maximum of the first burst, and |^^^|max is the maximum 
amplitude, including the contributions from the electromagnetic fields. Note that we speak of the 
inner core, where the matter falls subsonically, which corresponds to the unshocked region after 
core bounce. 



Model 




Pmaxb 


Mi.cb 


T/|VF|final 


-Em/ final 


At 


\h^^\ 

1 1 max 




(ms) 


(10^"^ g cm-3) 




(%) 


(%) 




(ms) 


(10-20) 


SSL7 


214.6 


2.97 


0.61 


1.0 X 10-2 


3.5 X 10" 


-8 


0.75 


3.5 X 10-^ 


MSL4 


242.8 


2.61 


0.83 


8.7 


4.3 X 10- 


-5 


0.74 


1.66 


MSL2 


242.9 


2.66 


0.83 


8.6 


4.2 X 10- 


-3 


0.73 


1.66 


MSLl 


242.7 


2.65 


0.83 


8.6 


3.9 X 10" 


-2 


0.73 


1.65 


MSS4 


241.8 


1.72 


0.87 


9.0 


2.9 X 10- 


-4 


0.58 


1.78 


MSS2 


242.0 


1.71 


0.87 


9.0 


2.9 X 10- 


-2 


0.58 


1.80 


MSSl 


242.1 


2.10 


0.87 


9.1 


2.4 X 10- 


-1 


0.58 


1.87 


MCS4 


245.1 


1.45 


0.91 


9.9 


3.0 X 10" 


-4 


0.53 


2.54 


MCS2 


244.1 


1.55 


0.91 


10.0 


2.7 X 10- 


-2 


0.53 


2.58 


MCSl 


244.2 


1.88 


0.91 


10.1 


1.8 X 10- 


-1 


0.53 


2.63 


MCS0.3 


245.3 


2.25 


0.91 


10.0 


5.8 X 10- 


-1 


0.53 


2.34 


RCSl 


302.2 


0.17 


1.11 


12.2 


3.1 X 10- 


-1 


2.03 


0.90 


RCS0.3 


320.5 


0.84 


1.14 


12.2 


8.0 X 10- 


-1 


1.76 


1.14 


MSL4-LS 215.4 


3.20 


0.69 


8.2 


3.4 X 10- 


-5 


0.63 


1.62 



equation of states commonly used in the supernova simulations. Since Lattimer-Swesty 
EOS is softer than the relativistic EOS, the mass of the inner core for the Lattimer-Swesty 
EOS becomes small, which reduces the mass quadrupole moments at core bounce. On the 
other hand, the softness allows the core to contract deeply, which makes the central density 
larger at core bounce. By the competition of these factors, the maximum amplitude remains 
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TABLE III: Effects of the strong magnetic fields on the maximum amphtudes. Note that ah the 
values of the amplitudes of the gravitational wave are given in unit of 10"^'^ in this table. 



Model 


"■quad "-jxE 








I^JJfi/^quidK^) 


MSLl 


- 1.65 2.40 X 10" 


-2 


-8.56 X 10- 




1.5 


MSSl 


- 2.16 2.95 X 10" 


-1 


-1.04 X 10" 


-3 


13.8 


MCSl 


- 2.69 6.46 X 10" 


-2 


-1.01 X 10" 


-3 


2.4 


MCS0.3 


- 2.55 2.13 X 10" 


^1 


-4.63 X 10" 


-3 


8.3 


RCSl 


- 1.01 1.11 X 10" 


-1 


-2.45 X 10" 


-3 


10.9 


RCS0.3 


- 1.53 3.83 X 10" 


-1 


-1.92 X 10" 


-3 


25.0 



almost the same between the two realistic EOS's. We note the difference that the typical 
frequencies of the gravitational waves for the Lattimer-Swesty EOS become higher due to 
the softness of the EOS. These effects of the EOS's on the gravitational waves are found to 
be independent of those of the magnetic fields, since the strong magnetic fields change the 
maximum amplitudes and do not change the typical frequencies of the gravitational waves. 

(3) For the weakly magnetized models, whose initial -Em/|iy|'s are smaller than 0.1%, the 
gravitational wave amplitude and the waveform are determined by rotation as in the pure 
rotation case j^^. The waveforms are categorized into the conventional criteria ^. The 
maximum amplitudes are within the detection limits for the detectors of TAMA and first 
LIGO for a supernova at a distance of 10 kpc. From the model based on the recent stellar 
evolution calculation, it seems difficult to detect the gravitational wave unless a supernova 
occurs very close to us. Type III does not occur if the two realistic EOS's are employed. 

(4) The signs of the values of the second peaks are found to be negative for the strongly 
differentially rotating models with a cylindrical rotation law, and positive for the other 
models. Due to the j x B part of the electromagnetic contributions, the positive values 
of the second peak become more positive, while the negative values become less negative. 
However the absolute values of the second peaks are within the detection limit of the first 
LIGO for a supernova at a distance of 10 kpc. It means that it will give us information 
about the angular momentum distribution of massive evolved star if a supernova occurs at 
our galactic center. On the other hand, it seems difficult to extract the information about 
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TABLE IV: Some characteristic quantities for the waveform analysis. The names of the models 
whose initial rotation law is cylindrical with strong differential rotation are written in bold letters. 
T^gj. and T^J^ the first and second oscillation period of the inner core, respectively, /i^^onj is the 
amplitude of gravitational wave at the second peak. 



Model 


rpl 

^ osc 

[ms] 


^ OSC 

[ms] 


second 

[10-20] 


MSL4 


1.9 


1, 


.3 


0.93 


MSL2 


1.9 


1, 


.1 


0.93 


MSLl 


2.0 


1, 


.2 


1.02 


MSS4 


2.4 


2, 


.8 


0.87 


MSS2 


2.5 


2, 


.7 


0.88 


MSSl 


2.5 


2, 


.8 


1.03 


MCS4 


3.1 


2 


.7 


- 0.62 


MCS2 


3.1 


2 


.7 


- 0.61 


MCSl 


3.1 


2, 


.7 


- 0.63 


MCS0.3 


3.2 


2 


.7 


- 0.57 


RCSl 


11.2 


9 


.1 


- 0.47 


RCS0.3 


11.4 


9 


.2 


- 0.53 


MSL4-LS 


1.5 


1.1 


1.00 



the magnetic field only by the gravitational wave from a magnetorotational collapse of a 
supernova. By the lack of the significant features by the magnetic fields in the gravitational 
signals in addition to those by rotation, a weak upper limit of the initial magnetic field 
strength of ~ 10^^ G in the central core prior to gravitational collapse may be put. 

Finally, we state some discussions based on the results in this study. In all the models, 
the final rotation rates are within the critical value where Maclourion spheroids become 
secularly unstable against triaxial perturbations (T/|iy|finai < 13.75%). Therefore, the axial 
symmetry assumed in this study may be justified. It is also noted that Rampp et al. 0, 
who performed the three-dimensional (3D) rotational core-collapse simulations, suggest that 
the features of the gravitational signals do not differ significantly from those by the two- 

17 



dimensional simulations. However, it is necessary to perform the 3D simulations for the more 
realistic estimation of the gravitational wave signals. As for the magnetic field contributions, 
the entanglement of the field lines in the azimuthal direction may suppress the gravitational 
radiation for a given rotation rate. Furthermore, it is noted that the non-axisymmetric 
motions of the matter produce the cross modes of the gravitational waves, which vanish in 



case of axial symmetry 
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44| . From the recent results by Fryer & Warren (2003) |4J| , who 



did 3D rotational core-collapse simulations, the maximum amplitude of the cross modes was 
presented to be ~ 10^^^ at several 100 Hz for a supernova in our galactic center. Although 
this value is smaller than that of the plus modes, the value seems to be within the detection 
limits of the gravitational wave detectors in the next generation, such as advanced LIGO 
and LCGT. If the characteristics of the waveforms between the plus and cross modes are 
detected, we think that it may be a good tool to extract information not only about rotation 
but also about magnetic fields in the evolved massive star. In order to investigate this, we 
are currently preparing for the 3D MHD simulations (H. Sawai et al. 2004 in preparation). 

Recent observations imp lv that gamma-ray bursts (GRBs) are associated with core- 
collapse supernovae j3, lij- Although the failed supernova or the so-called collapsar model 
are supposed to be probable 47|, |48[, it is still controversial what powers GRBs. As the 
energy source for fireballs, the energy deposition by neutrinos emitted from the accretion 
disk or by the MHD process which extracts the angular momentum of the rotating black 
hole are considered to be important. MacFadyen et al. [4^ reported that matter with the 
mass of ~ 0.1 — 5M0 falls into the central black hole during minutes to hours, which suggests 
the duration for the activity of GRBs. It is noted that no magnetic fields are considered in 
their calculations. When we boldly extrapolate the central magnetic fields in our strongest 
magnetized models, which become as high as ~ 10^'' G in the central region, to the model 
of GRBs, it may not seem so unnatural to imagine the infall of the strongly magnetized ma- 
terial into the black hole. If the magnetized matter falls into the black hole anisotropically, 
we expect the emission of the gravitational waves. In order to calculate the amplitudes, the 
formula, which we derived in this paper, will be useful. Here it should be noted that the fully 
general relativistic magnetohydrodynamic simulations with the multi-dimensional neutrino 
transport are required for the reliable estimation to determine the masses of the black hole 
and the accretion disk. Therefore further advancements in the numerical simulations are 
necessary for the quantitative discussions. 
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Recently, one group claims the detection of linear polarization of ~ 80 % in a GRB 
Isoj ]. however, the other group claims that the polarization is less than ~ 4 % 51 1. The 
discrepancy may come from the difference in the way of data analysis. If the former is 
correct, it may seem consistent with our results, which predict rather coherent magnetic 
fields in the vicinity of the central core. 

As stated earlier, Miiller et al. jl^ recently reported that the maximum amplitude is not 
obtained at the core bounce but at the later phases when the anisotropic neutrino radiation 



becomes active (see, also, 
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43|). If the degree of the anisotropy of neutrino radiation 



could be observed from the gravitational radiation, it will help to understand the explosion 
mechanism itself. This is because anisotropic neutrino radiation induced by rotation is likely 



to p roduce an asymmetric explosion 



to p roauce an a 

y, 55, Q 



33, 



44l . l52j | as suggested by observation of SN 1987 A 



Unfortunately, this is beyond our scope now to estimate the effects of 
anisotropic neutrino radiation on the gravitational signals by our crude treatment of the 
neutrino transport. For the purpose, we are now developing a two-dimensional neutrino 
transport code, which is indispensable for the estimation of the gravitational wave from the 
anisotropic neutrino radiation. 
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FIG. 1: Waveforms for models MSLl (left panel), RCSl (right panel). In the left panel, "quad", 
"j X B ", and "rhom" (pm) represents the contribution from the mass quadrupole moment, jxB part, 
and the time derivatives of the magnetic energy density (pm part), respectively. The amplitudes 
are artificially multiplied by 10^ for j x B part and 5 x 10'^ for pm part to make their waveforms 
clear. From the panel, j x i? part is found to dominate over pm part. In the right panel, "tot" 
indicates the amplitude including the total contributions. From the panel, it is found that the 
j X B part changes at the opposite phase of the mass quadrupole moment. Note that the source 
is assumed to be located at the distance of 10 kpc. 
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FIG. 2: Waveforms for models MSL4 (left panel), MCS4 (right panel). Note that the source is 
assumed to be located at the distance of 10 kpc. 
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FIG. 3: Relation between the central density and the effective adiabatic index 7 at near core 
bounce. Relativistic EOS, LS EOS in the figure represents the relation taken from the values from 
the models MSL4 and MSL4-LS, respectively. The swirls at the density of ~ 10^^'^ g cm"^ are due 
to the reexpansion of the core after core bounce. 
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FIG. 4: Waveforms (left panel) and the time evolution of the central density (right panel) for 
models MSL4 (left panel), MSL4-LS (right panel). All times are relative to the time of bounce {%) 
in this figure. Note that the source is assumed to be located at the distance of 10 kpc. 
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FIG. 5: Waveform for model SSL7, which is based on the recent stellar evolution calculation. The 
amplitudes are about three orders of magnitudes smaller than the other models. Note that the 
source is assumed to be located at the distance of 10 kpc. 
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FIG. 6: Detection limits of TAMA [Sfl], first L 
Cryogenic Gravitational wave Telescope (LCGT) 
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, and Large-scale 



with the expected amplitudes from numerical 



simulations. The open squares represent the maximum amplitudes for all the models, except for 
models MSL4 (closed triangle) and MSL4-LS (closed circle). On the other hand, the pluses and the 
closed squares represent the amplitudes of second peaks for models with strong differential rotation 
with a cylindrical rotation law and for the other models, respectively. Open triangle represents the 
maximum amplitude for model SSL7. We estimate the characteristic frequencies by the inverse of 
duration periods of the corresponding peaks. Note that the source is assumed to be located at the 
distance of 10 kpc. 
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